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The mass spectrum of the primordial black holes formed by density perturbation in the radiation- 
dominated era of the Randall-Sundrum type-2 cosmology is given. The spectrum coincides with 
standard four-dimensional one on large scales but the deviation is apparent on smaller scales. The 
mass spectrum is initially softer than standard four-dimensional one, while after accretion during 
the earliest era it becomes harder than that. We also show expected extragalactic diffuse photon 
Ch , background spectra varying the initial perturbation power-law power spectrum and give constraints 

^ ■ on the blue spectra and/or the reheating temperature. The most recent observations on the small 

scale density perturbation from WMAP, SDSS and Lyman-a Forest are used. What we get are 
, interpreted as constraints on the smaller-scale inflation on the brane connected to the larger one 

at the scale of Lyman-a Forest. If we set the bulk curvature radius to be 0.1 mm and assume 
' the reheating temperature is higher than 10® GeV, the scalar spectral index from the smaller scale 

^ ' inflation is constrained to be n < 1.3. Typically, the constraints are tighter than the standard 

four- dimensional one, which is also revised by us using the most recent observations. 

o 
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I. INTRODUCTION 

Q . It is well believed that the life of heavy stars ends with luminous explosion coUapsing into a black hole. The 
}—{ ' birth of such massive black holes are daily phenomena in that we can observe supernovae at the present time. On 
c/j . the other hand, the mechanisms to produce holes smaller than stellar mass in the early Universe had been invented 
' several decades ago 0, Q- Those very small holes, called primordial black holes, are thought to be formed from 
density perturbation in the form of radiation, and their weight is comparable to the horizon mass in the radiation 
dominated era The primordial holes can exist everywhere in the Universe. Those in the dense regions might 

have clustered, being one of the candidates for the dark matter. The significance given to the primordial holes was 
not only their existence itself, but also the possibility of the observation of the Hawking radiation . The quantum 
■ ■ ' theory on the black hole background predicts that black holes radiate particles obeying almost black-body spectrum 
0- It also predicts that the stellar mass holes cannot produce particle fields significantly, whereas primordial black 
holes are likely to emit various heavy particles, some of which may reach to the earth as cosmic rays. One of the 
earliest work in this respect was done in |3| , which considered the possible contribution of PBHs to the extragalactic 
gamma-ray background. Constraints on the PBH abundance from cosmic rays have been extensively investigated 
Hi ^1 ■ Primordial black holes serving as the sources of cosmic rays or relic particles were also studied in the 
literature [Hill. A turning point for PBHs comes with the discovery of the cosmological inflation [11 m. They 
predicts almost Gaussian, super-horizon scale density perturbations which are thought not only to serve as seeds of 
the structure formation in late time, but an origin of primordial black holes. Then, in addition to the large scale 
structure like the cosmic microwave background (CMB) or the spatial distribution of the galaxies, PBHs became a 
probe for rather small scale structure in the inflationary early Universe [ll 111 111 111 l20l |2U I2I I2I l25l| . 

After some decades, the study of cosmology is shifting toward higher-dimensional models called braneworld. The 
idea of a membrane as our Universe has been inspired by the string theory. A toy model for the fundamental 
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theory suggested by Randall and Sundrum |26l l27l | is one of those which successfully matches with standard big- 
bang cosmology in the late time, but deviates in the earliest epoch of the radiation dominated-era |2^ Isol l3ll | 
(see, for example, 32] for a review). In the context of such modern cosmology, we should rediscover the utility of 
primordial black hole. Particularly, PBHs in the Randall-Sundrum (RS) cosmology show very different features from 
those in standard four-dimensional scenario 33, 34, 35. 36] . They covers all the importance, i.e., formation dynamics, 
accretion process, and the Hawking radiation. Primordial black hole is an exceedingly good probe for the brane early 
Universe. Possibilities for constraining the braneworld geometry, such as the bulk curvature radius, have been shown 
by investigating the cosmic-ray spectra 37, 38]. 

In this paper, we give the connection between the curvature perturbation spectrum generated on the brane after 
brany or bulky inflation (for example, j39l l40l l4l|') and the resulting mass spectrum of the primordial black holes. 
In addition, we constrain the blue spectral index of the curvature perturbation using the observational data of the 
extragalactic diffuse photon background. The cosmological parameters we use are based on the concordance model 
containing the cosmological constant and cold dark matter (ACDM), which has been recently confirmed by the 
observation of cosmic microwave background (CMB) by Wilkinson Microwave Anisotropy Probe (WMAP) p3. l43l]. 
For cosmological parameters, we use First Year WMAP and Two Degree Field system (2dF) data for the density, 
Hubble parameter, age of the Universe, and so on. On the other hand, for smaller-scale matter power spectrum, 
we use the result of jij], which combined WMAP data with Lyman-a forest (LyaF) and Sloan Digital Sky Survey 
(SDSS). With these recent observational data, we also reexamine the constraint on the perturbation spectrum in four 
dimensions [TqI ]. 

This paper is organized as follows. In Sec. Ill Al we review the basic framework to compute mass spectrum of 
gravitationally bounded objects like primordial black holes from density perturbations by the aid the Press-Schechter 
theory. Then after describing the Randall-Sundrum braneworld in Sec. IIIBI we show the mass spectrum of primor- 
dial black holes in the RS cosmology. Finally in Sec. IIIII we give constraints on brane inflation via the curvature 
perturbation spectrum and the reheating temperature. 



II. MASS SPECTRUM OF PBHS IN THE RS BRANEWORLD 
A. Press-Schechter prescription for inflationary density perturbation 

In this section we summarize how to derive the mass spectrum of PBHs in the framework described by Press and 
Schechter |4^]. In the four-dimensional theory there has been much effort to obtain the mass spectrum of PBHs. Since 
we are now focusing on possible differences emerging from different dimensionality, the mass spectrum is described in a 
manner such that the differences will be clarified. We assume the bulk spacetime has five dimensions while the matter 
content is confined in a four-dimensional brane. Although we will only consider the five-dimensional Randall-Sundrum 
model, the following arguments will not strongly depend on the detail of the model. 

First we start from the comoving number density of density fluctuation which arose in the inflationary early Universe. 
For the moment we do not care about the mechanism of the inflation in the braneworld, which would be eventually 
included in the power spectrum of scalar perturbation. Now we focus on the quantities evaluated at some fixed time 
in the radiation-dominated era. The fiuctuation is, after smoothed by some suitable window function, at every point 
spherical with comoving radius R, which is the scale of the window function. On comoving hypersurfaces, density 
fluctuation at x is characterized by 5(x) = (/5(x) — where p, without explicit spatial dependence, represents the 
background energy density. S is averaged by a window function W to give a smoothed quantity 

^fi(x) = ^J VK(]x' - x]/i?)<5(x')d3x', (2.1) 

where V = J VF(]x' — x]/i?)d'^x' is the comoving volume contained in the window function. The spectrum of the 
smoothed density contrast is given through its Fourier transform, Sr = (27r)'^/^ / e^^^'^dnd^x, by taking the ensemble 
average as Vs — [fc'^/(27r^)]()(5_R)^). Vs, assumed to be independent of the direction of k, gives the expression for the 
variance 

4 = ('5|(x))=_| W'{kR)Vs{k)-, (2.2) 

where W is the volume-normalized Fourier transform of W. Below we will work throughout with those smoothed 
quantities Sr, gr, and Vs, so the indication of smoothing, R, will be dropped hereafter. We introduce M = pV , 
the mass contained in the window function. The comoving number density of smoothed fluctuations whose comoving 
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radius and amplitude range within R E [R,R + dR] and S E [S,S + dS] , respectively, is written as 

n{R, 6)dRdS = ^F{R, S)dRdS, (2.3) 
where FdRdS is the fraction of the space. We assume the Gaussian probability distribution of the density fluctuation 

where the prefactor 2 is of the Press-Schechter prescription. It gives the fraction F as 



F{R,S) = 



dP 



dR 



2 Ict'I / S 
TT exp 



2 i n 2 



X2 



6^ 

— -1 



(2.5) 



where a prime denotes derivative over the scale of smoothing, R. 

Next we describe the spectrum of density perturbation. Usual inflationary models predict that curvature pertur- 
bation on comoving hypersurfaces, 7?,c, is produced during accelerating expansion, and its power spectrum V-ji^ tends 
to obey power-law. Cosmological perturbation theory tells us how such a perturbation evolves and what we get as 
observables. Now we want to know the spectrum of density perturbation, which form PBHs when they cross inward 
the cosmological horizon. In the contexts of higher-dimensional cosmology, there are certain points to be taken se- 
riously. We first notice that the presence of the extra dimensions may highly affect the gravitational dynamics of 
scalar perturbation on a brane. The effects can be categorized into two classes; one is locally determined by matter 
on the brane while the other is the non-local responses of the bulk geometry to the perturbation on the brane. The 
local effect coming from non-standard form of matter source is reflected in the background cosmological expansion, 
which of the Randall-Sundrum cosmology we review in the next subsection. To see how the latter affects, let us take 
a look at the evolution of scalar perturbation in the RS cosmology. In the super-horizon regime, above introduced 
T^c coincides with another variable C^^ the gauge-invariant curvature perturbation on uniform-density hypersurfaces 
|4g . Four-dimensional inflationary theories predict that these values are determined by the inflaton expectation- value 
fluctuation and are almost constant outside the cosmological horizon On f^e other hand, in the five-dimensional 
Randall-Sundrum cosmology, a non-adiabatic Kaluza-Klein mode contributes to ^ and generally it varies with time 
|32, 48, 49, 50J, which can modify the resulting matter spectrum. In addition, the initial value measured at the exit of 
horizon during inflation is also modified by the higher-dimensional effect [5ll|52,|53- We should include such non-local 
effects from the bulk space for the complete arguments, however, to solve perturbation in the extra dimensions is 
difficult and, in particular, discussions for determination of initial scalar perturbation |5ll l52l Is^ seems not fully 
satisfactory. Therefore, in this paper we put apart such non-local effects but focus on the local background change. 
In the language used in for the RS cosmology, what we do is to set those contributions from anisotropic stress, 
dark radiation, and isocurvature originating from higher-dimensional nature to be zero; nonetheless, we are still able 
to capture one of the essential features of higher-dimensional cosmology. The constancy of curvature perturbation 
and its power spectrum in the super-horizon regime is recovered in this framework. A relevant brane inflation model 
is, for instance, what was treated in 39]. 

Now let us derive and investigate the density perturbation spectrum. At sufficiently late times, i.e., four-dimensional 
regime, the comoving curvature perturbation 7?.c effectively takes a role as gravitational potential and induces observ- 
able density perturbation, whose spectrum Vs on large scales is given in the following manner |47| 

2u;\V fc 



where k = 1/R is the comoving wave length, w = p/p represents the equation of state of the era considered, and 
the transfer function T{k) expresses the effect of deformation to matter perturbation. H is the Hubble and a is the 
scale factor normalized to be unity at the present. The time when the above spectrum is evaluated can be arbitrarily 
chosen if only the Universe can be regarded as four dimensional then. Here we set the time by the condition that the 
Hubble radius becomes equal to some length a/fco- Then the variance is expressed as 

r°° dk 

/ W\kR)rs{k,t{aH ^ ko))^. (2.7) 

Through the above relations, the observed density perturbation spectrum J^^"^^^ can be utilized to constrain the 
curvature perturbation. Parameterizing as 

"a-l 



"'pivot 
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a combined result from First Year WMAP, SDSS galaxies (SDSS-gal), and Lya Forest claims Ag = 2.45 x 10~^ 
and Us = 0.971 at the pivot scale fcpivot — 0.05 Mpc~^ This parameterization is valid down to LyaF scale. 

Its corresponding wave number is a rather useful choice as the reference scale for our purpose, i.e., we now take 
ko = 0.9h Mpc~^. This length is in the radiation-domination regime but is much larger than the scale at which 
the cosmological evolution is possibly modified. Since the scales we wish to consider are much smaller than those 
appeared in the observations mentioned above, we now assume a single power-law spectrum below the LyaF scale 
reparameterized as 

'PnAk) = A -) , (2.9) 



ko 

where the amplitude is A = As{ko / kpi^ot)"" ~^ — 2.28 x 10~^, but now that the spectral index n is regarded as a free 
parameter. In other words, what is assumed here to realize such a broken power-law is some generation mechanism 
like the double inflation model where the connecting scale between the first and second inflations is placed at the 
Lya Forest scale. From this point of view, what we will consider can be comprehensively interpreted as limits on the 
second inflation. Although not all the brane inflation models prefer such a fluctuation spectrum, we should take it 
into account that an enhancement of small scale density perturbation is observationally inferred by the First Year 
WMAP result jHlil. 

Substituting the spectrum into Eq. (|2.6() and setting w = 1/3 give the explicit form for the variance as 



>VM(fcoi?)"^^+"^ = (T^(fcoi^)"^^+"^ (2.10) 



CSO 

W''iR)= I W^{x)T^{x/R)x^+''dx, (2.11) 



where we defined W as 





and (Tq = (4/9)'^W^A. If we choose the Gaussian filter W = ^ and, for simplicity, assume T — 1, then one 
gets = (l/2)r[(3 -I- n)/2]. The volume inside the Gaussian filter is V ^ {VOjiRf. 

Then we consider gravitational collapse occurring at the horizon entry. We denote the corresponding time as th- We 
introduce "linear evolution" factor D by which the amplitude and root-mean-square of density perturbation crossing 





inward the horizon are related to the values at the reference scale fcn as^ 



^'^Uz^-i/fj. (2.12) 



The following discussion does not depend on the detailed form of _D, which we will give in the next subsection. The 
collapsed black hole mass, Afbh, is determined by the amplitude (5h and the comoving scale R. Then the comoving 
mass spectrum of primordial black holes formed by density perturbation is now given as 



dnbh '■^'^ 



n(i?,5)J(Mbh,<5h)d4, (2.13) 



dMbh Js________ 

where the domain of integration is determined for each collapse model and 

J(Mbh,<5iO^ .f/f'^]. (2.14) 

is the Jacobian matrix. The integration is performed with fixed Afbh- Substituting Eq. H2.3|l and defining da: = d(5h/CTh, 
the above mass spectrum becomes 

dnbb [2 r- J (2.15) 



dMbh V TT ./^,„,„ VD a 



^ The time sequence t{aH = ko) > tji is understood by the negative power of D. 



5 



In usual cases, the scale R determines the time of horizon entry ih- This deduces the functional form 5h to be 

(5h = D{R)-^d. (2.16) 

The collapsed mass Afbh may be a function of both the comoving scale R and the amplitude (5h. In ordinary parame- 
terization, the form of Mbh can be expressed as 

Mbh = /((5h)A4(i?), (2.17) 

where Mh is the mass contained in the Hubble volume, the horizon mass, now being evaluated at the time of horizon 
entry, and / is the fraction. Then the Jacobian is found to be J = D/M^^ and the formula for the mass spectrum 
reduces to 



dribh /2 /■^"'"' 1 



dMbh V TT J^^^^^^ VM(^^^ 



l\dx. (2.18) 



Gravitational collapse taking place at the horizon-entrance is an issue even in four dimensions, for instance, some 
critical phenomenon has been found |54l | . The process may also be affected by dimensionality. One simple prescription 
is to regard the emerging black hole mass almost the same as the horizon mass at the horizon entrance, i.e., / = 
const ~ 1. Actually Carr argued that / — w^^^ where w = p/p Q. For this reason, we here assume that / is 
almost constant and will concentrate on the relation between the comoving scale and horizon mass. Four dimensional 
constraints from the critical collapse model have been considered elsewhere IH^ . 

We estimate the PBH mass spectrum a bit more assuming constant /. Another assumption that both Afbh and a 
are power functions of R gives \a' \ / aV M^^^^ ~ 1/yAfbh, so that one obtains 

L_ r^^^^ e~-"/^\x^ - l|dx = -i— I(Afbh). (2.19) 



dAfbh V^Afbhy:.„„„ V^Az/bh 
The integration is performed to give 

{-^(*^min) -^(-^max) 

-/^(a;max) - I{Xmin) X^a^ < 1 , (2.20) 

2/(1) - /(xmax) - lixmin) Otherwise 

where I{x) = xe~^ 1"^ . Now we observe that the denominator yAfbh gives the "scale- invariant" part of the mass 
spectrum that does not depend on the density perturbation spectrum. On the other hand, the integral X determines 
its scale-dependent tilt through In general (5h,min ~ <5h,max ~ 0(1). If we adopt power-law dh oc A:" normalized on 
the large scale k^^ as fih ^ Vj4 ^ 0(10^^), then we deduce that the PBH abundance never possesses finite values 
unless (Jh at least grows up to <5h,min ~ 0(1). Hence, hereafter we consider such blue spectrum cases with positive a, 
where the explicit form of a expressed with spectral index n will be found later. In such a situation, there would be 
a threshold scale i?th at which (Jh ^ (5h,min- T- is estimated as X ~ ^^mine^^™'"^^ above the threshold, while X ~ Xmax 
below it. Therefore the approximated expression for the tilt is given as 

,exp 



ul\R) X I '^"'-'"-^ L ^<^R)\ ^ - . (2.21) 

l.'5h,max R ^ ^th 

We can also evaluate the mass density of the PBHs from the mass spectrum which has been often used in the 
literature. phh{> Af ), the mass density of black holes whose mass is greater than Af , can be evaluated as 

Pbh(> M) = r A/bh^dAfbh - r §^^dAh^, (2.22) 
Jm dAfbh Jm y(Mhh) 



which further leads to 



^-.An ^ jx^in(Af)-ie-™W/2 Af>Afth 
^-(>"^)-nM)M-n.ax(A/) A.<A/,b' ^'-''^ 



where Afth is the mass scale associated with i?th- This phh is dominated by the holes whose mass is M because 
the mass spectrum is a rapidly decreasing function with increasing Afbh- It is useful to define the mass fraction 
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of black holes relative to the background. The fraction of black holes whose mass is (larger than) M is given as 
P{M) = phh{> M)/ P- This is a time-dependent quantity. For convenience, we use /3i, which is the ratio to the 
background energy density at the time of the formation of holes with mass M . Since such a radiation energy density 
is given by pi = M /V , the initial fraction is obtained as 



A(M) 



(Jh (^^ ) exp 



M > Ma, 
M < Mth 



(2.24) 



We note that with the above assumptions this (3i is equivalent to a, defined as the ratio of mass density of holes relative 
to radiation energy density at the time of formation, which has been also used in the past works [itL Is^ fssl IH^ . 

In the following our task is to obtain expressions for V, Mbh, and dh according to each cosmological model. We 
then will be able to constrain inflationary models through the perturbation spectrum by constraints on the PBH 
abundance. 



B. Randall-Sundrum braneworld cosmology 



In this section we describe the background cosmology and the evolution of density perturbation in the so-called 
Randall-Sundrum type-2 braneworld pa. l27l |. This theory introduces small-scale modification whereas the late-time 
behavior converges to the purely four dimensional one psL I29I Isol Isif . Hereafter we basically work with the units 
such that all the four-dimensional Planck units are unity: M4 ^ T4 = £4^ — — I, although sometimes they appear 
in expressions. 

The RS2 braneworld is composed by a four dimensional timelike brane, which has positive tension, embedded in 
five dimensional anti de Sitter (AdS) bulk. Standard model matter is confined within the brane. The cosmological 
expansion on the brane is determined by the modified Friedmann equation 



, Svr 



3Mi 



A4 K 

T 



(2.25) 



where A > is the tension of the brane, ps the energy density of the so-called dark radiation, and K the topology of 
the three dimensional space. The effective Planck scale M4 and effective cosmological constant on the brane A4 are, 
respectively, given as 

where M5 is the fundamental mass scale in five dimensions and A5 < is the cosmological constant in the bulk. 
It is useful to define the scale of the theory by the curvature radius of bulk AdS £ = y'iMf/{A^n\K^ . Setting A4 
sufficiently small, the relation between the geometrical scale and the brane tension is obtained as A = 3M4/(47r^^). 
Then we get M| = £M|. In the following we omit A4 and K , which have been confirmed to be sufficiently small by 
WMAP, and also drop p£ js^ • Then it is obtained that the cosmological expansion in the radiation dominated era 
is separated into two parts, namely in the earlier one the p^/A term dominates ("p-square era") whereas in the latter 
the p term does. We often denote the two separated phases simply by "p^" and "p^"j respectively. For later use, we 
define the effective energy density and pressure 

Peff = P (1 + 1^) ' (2-27) 
Pcff-P+^(2p + p), (2.28) 
which give the modified Friedmann equation Eq. H2.25|l . We can represent the effective equation of state as 

-off-^^jf (2.29) 

Note that in our formulation, all the expressions and corresponding results associated with the p^ case completely 
agree with the standard four-dimensional cosmology since the £ ^ limit recovers the 4D cosmology. 
The solution through the whole radiation-dominated era is 



^1/4(^ + ^^1/4 
tl 



a{t) = a,q -TJT^^ (2.30) 
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where tc = £/2 and the normaHzation is given at the time of matter-radiation equahty toq. The Hubble parameter is 
given as 



Hit) 



2t + tc 



(2.31) 



4t(t + U) 

Other background quantities, radiation energy density p and the horizon mass = AirpKiH^), are calculated as 



p{t) = 



AtiP 
1 



^/mi^ + 1-1 



327ri(t + ic)' 



{2t + tc 



(2.32) 
(2.33) 



It is useful to note the asymptotic behavior in respective regime: Mh oc t^ for and Mh oc t for p^. Knowing the 
background cosmological evolution, the horizon-crossing time for a density fluctuation with radius R = 1/k is 
determined by a condition 



Ra{ti,) = H-\ti,) 



(2.34) 



Now let us derive the connection between the quantities at some fixed time (i?, d) and those at the time of gravitational 
collapse (Mh, (5h). The comoving radius R determines the horizon crossing time th from the condition 



R 



tii' K^\t^+tcr/' 



2th + tc 

while the horizon mass at ih is evaluated as Afh(th)- Eliminating t^ from these expressions gives 

1/2 

I / Mj^'-M^" \ 

II 71 T I 

koR'. 



(2.35) 




(2.36) 



where Mh,o = M\i[t{aH = ko)] « 1.94 x 10^^ g is the horizon mass at our normalization scale and we defined 
Mc = £/4 « 3.37 X 10^^(£/0.1 mm) g. One can adopt the above relation to the primordial power-law perturbation 
spectrum to find 



Cr ~ CTg X < 



3/4 .,1/4 



-(3+n)/4 



-(3+n)/4 



(2.37) 



Next, we consider the evolution of scalar perturbations. We focus on super-horizon regime and, for the reasons 
described in the previous subsection, omit anisotropic stress, dark radiation including its p erturbation, and isocurva- 
ture. Then it follows that the gauge-invariant curvature perturbation C conserves 'is, 49, "SO^, which leads to, in the 
longitudinal gauge metric perturbation $ (now equals to —'^) is also constant (49.j : 



^ _ 3 + 3wcff ^ 
5 + 3woff 



(2.38) 



On large scales C coincides with TZc ■ From the remaining perturbed Einstein equation, the relation between $ and 6p 
to the linear order is found as 



An 
Ml 



1 



(2.39) 



where a nabla denotes vector derivatives over the comoving spatial coordinates. From the above relations, we even- 
tually find the time evolution of <5, or equivalently that of its Fourier transform 6, as 



jL (Jl\^7? 

15 \ aH) 



4 / fe 
. 9 



(2.40) 
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The latter was used to give Eq. p. 6(1 . The evolution factor D is then obtained as 



a 



5h 



Finally we get 



(To X 




(l-«)/4 



(2.41) 



(2.42) 



Throughout the above arguments, we have implicitly assumed that the beginning of the radiation-dominated era is 
before the transition. Otherwise we do not see the effect of the fifth dimension. The time of the beginning is regarded 
as the end of inflation, namely, reheating. The reheating time irh is expected to satisfy irh < • The scale associated 
with reheating will be discussed in the next subsection together with the PBH mass spectrum. 



C. The mass spectrum of PBH 



Now we have obtained all the relations connecting the initial PBH mass Mbh and comoving scales such as R 
through the parameters specifying cosmology, i.e., in this case £. Here the word initial means the spectrum is given 
as a function of the mass soon after gravitational collapse. In this sense, we denote it as Alhh.i- Again we set / to be 
constant, then the mass spectrum is written down as 



dn 



bh 



dA/, 



bh.i 



2 3 + n PoqMcV^ 



-3/8, ,-17/8 



(1 



i) 



o-h(Mbh,i) 



-5/2 



bh,i 



p'(Mbh,i < /Me) 
P'(Mbh,i>/Me) 



(2.43) 



where 




(n-l)/4 



p'(Mbh,i < /Mc 



(2.44) 



We defined poq = p(ioq) and M^q = pcqV{Rcq)/il + z^qf ~ Mh(toq) ^ 10^^ g. For (5h,min, we refer to 0, in which 
the collapse threshold on the background was found to be (5h,min ~ 0.1 rather than 1/3 of the ordinary radiation 
dominated era 7|. The mass spectrum coincides with the four-dimensional one at larger scales, while on smaller scales 
it is deformed to be soft. In the RS cosmology, we further incorporate an additional effect on the evolution of PBHs, 
which is not the mass loss via Hawking radiation. On the contrary to the ordinary 4D case, radiation fluid accretes 
onto PBHs efficiently in the p-square phase due to the slow expansion rate, which causes power-law growth of the 
mass Afbh(0 i^''^/'^ j34|, |35| , where < F < 1 is the efficiency of accretion (for detailed treatment of accretion and 
its consequences, see |53). Only smaller holes formed in the p-square phase suffer accretion soon after their formation 
till the end of the phase. Since PBHs now being interested in are those which survive until at least the epoch of 
nucleosynthesis, we can ignore their mass loss during p-square phase. Then the mass at the end of accretion, which 
we denote primordial mass Mbh.p, is given as 



Mbh,p = I — 



2F/TI 



bh.i 



(16/A4)^/"A^(^^)/^ 



(2.45) 



Inverting this gives A^bh.i = (IG/ATc)"-^/^''"-^^ Af^/^''"-^^ 



mass Afbh.p is given by 
d^bh _ dAfbh,i dnbh 



bh,p 



dAfbh.p dA/bha 

+ n p^qMll^ 4,, 



The PBH mass spectrum with respect to the primordial 

(2.46) 



^£(Afbh,i) 



(1 + ^^cq) 



(Th(Af, 



bh,ij 




^3/2^,^-5/2 



bh,p 



bh,p 



p2(Afbh,p</A/c) 
pHMbh,p>/Me) 

(2.47) 
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where T = 9F/[8(7r — F)]. Mass spectra for £ = 0.1 mm and / = 1 are shown in Fig.^ The left panel shows n = 1.00 
case (scale- invariant), whereas the right is for the n = 1.60 case. In each panel, accretion efficiency for two 5D spectra 
is taken to be F = 1 and F = 0. 




Mbh.p [9] Mbh,p [g] 

FIG. 1: Left panel: n = 1.00 case, or the scale invariant part, of the mass spectrum. The spectral index for standard 4D is 
—5/2. 5D spectra coincide with it above the threshold fMc ~ 10^^ g, whereas below that, their indices deviate from —5/2: 
— 17/8 — —2.125 for F = and —2.65 for F = 1, respectively. _F = spectrum is also regarded as the initial mass spectrum 
before accretion. Values on the vertical axis are only expressing the relative scale. Right panel: n = 1.60 case. Heavier holes 
than A/th are exponentially suppressed. 



We note that the scale of reheating sets the small-scale cut-off in the spectrum. In our treatment, the horizon mass 
at the reheating Mh,rh = -Mh(irh) gives the minimum mass in the spectrum. The initial condition after inflation is 
comprehensively given via the reheating temperature T^h. In the p-square era the modified Friedmann equation gives 
cx p^/X oc g^T*£^, where we used the ordinary expression for radiation energy density p = (tt^ /30)gT* with g the 
particle degree of freedom. Hence horizon mass at the reheating time is 
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Mh,rh= {^^j g~''rX-^ (2.48) 



^2.37xl03M^ g, (2.49) 

VlOO/ VO-lmmy VlO^GeVy ^' ^ ' 

which corresponds to the smallest PBH mass Afbh.min = /-^^h,rh- If one wants the Universe to begin in the p^-era, 
then Mji^rh ^ is required. This claims the minimum reheating temperature for the existence of 5D PBHs 

/ e 

Trh > 10^ — GeV. (2.50) 

\0.1 mm J 

In the 4D case, (x p cx gT"^, so 



MKr.-\l^9-'^Xl' (2.51) 



III. CONSTRAINT ON THE PERTURBATION SPECTRUM FROM PBHS 



We are interested in giving constraints on the density perturbation using limits on the PBH abundance. In our case, 
it will lead to informations on inflation in the RS braneworld. We use the large-scale normalization erg given through 
As obtained from WMAP+SDSS-gal-|-LyaF [4J|, and assume that the power spectrum below the LyaF normalization 
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scale obeys power-law with index n. Although the index Ug has been observed nearly scale-invariant at the LyaF 
scale, we cannot assert that it also applies on small scales. Hence, the footing where we stand in this section is that 
we regard the spectral index as a free parameter. Then blue spectrum with n > 1 can produce PBHs on small scales. 
Even if it is not the case, there is little reason to believe the absence of small-scale large density perturbation, which 
could be produced by some inflation mechanisms like double inflation. 

Overproduced PBHs generally cause various problems in each epoch of the Universe, hence their abundance at each 
scale is required to be sufficiently small. Among known constraints on PBH, the most stringent ones are set by cosmic 
rays assumed to be emitted via the Hawking-radiation process. At each early epoch in the cosmic history, the PBH 
abundance is constrained so that energetic cosmic-rays from holes must not spoil the standard big-bang cosmology, 
for example, radiation entropy [Glj . nucleosynthesis |62l |. and so on. Also in the context of braneworld cosmology, 
interesting issues arise from cosmic rays; the five-dimensional nature should be imprinted in the cosmic-ray spectra 
emitted from PBHs smaller than the anti-de Sitter curvature scale £. In this paper, we focus on a directly observable 
cosmic-ray in the present epoch, the extragalactic diffuse photon background. The most dominant contribution to 
the spectrum comes from those holes currently evaporating. Because of the red-shift dilution, contribution from the 
short-lived ones is suppressed, whereas low temperature of long-lived ones makes their emission power low. The upper 
bound on £ has been set around ~ 0.1 mm from gravitational experiments _62i, j6^J. As we will see, a lower border 
exists around ^ 10^^^ m, below which PBHs whose lifetime is equal to the age of the Universe have larger radius than 
the fivc-dimcnsional scale. In this case, five-dimensional effects on the observable cosmic-ray spectra are expected to 
diminish. 

The strategy is summarized as follows, (i) First we should determine of which size black holes take responsibility 
for the cosmic-ray spectra under consideration. In the RS cosmology, the mass scale is affected by the size of the bulk 
curvature scale i. (ii) Then the abundance of such holes is constrained by the cosmic-ray observations, (iii) From it 
we obtain constraints on the density perturbation using the mass spectrum formula. To formulate cosmic-ray spectra, 
we describe Hawking evaporation of PBHs in the next subsection. 



A. Braneworld black hole and its evaporation 



In braneworld scenarios with some large compactification scale like Randall-Sundrum, a massive object localized on 
the brane within a typical scale should see the background space-time as an effectively fiat five-dimensional Minkowski 
space (for possible deviations and specific exotic features, see [63, [63, 1^ IzS Ell)- Therefore a 

sufficiently small black hole is no longer ordinary 4D one, but is approximately described as 5D solution. If there is 
no charge and is not rotating, the solution should be five-dimensional Schwarzschild. The 5D Schwarzschild metric is 
given by 





['-(t)1 




['-(t)1 


-1 

dr" 













r^dnl 



(3.1) 



where rg is the Schwarzschild radius and dfln denotes the line element on unit n-sphere. The four-dimensional 
geometry seen by brane matter is obtained as its projection, i.e., after replacing df^s with df22- In our framework 
the Schwarzschild radius in the five-dimensional spacetime is related to the black hole mass Mbh and five-dimensional 
fundamental mass scale M5 as 



\3tt Ml J \3Tr\i4 



V Ma 



(3.2) 



while in four dimensions rg = 2(Afbh/Af4)i?4. The Hawking temperature of a Schwarzschild black hole is given in 
terms of its radius as Th = [D — 3)/47rrg, where D is the number of dimensions 74] . Thus in five dimensions 



27rra 



_3 

32"7r \li 



-1/2 / j.,f \ -1/2 
Mi 



(3.3) 



while in four dimensions Th l/(47rrg) = (l/87r)(Mbh/M4) ^T4. For later use we define the effective Schwarzschild 
radius for massless fields 



?'S,cff = 



D - i\ 1/(^^3) / D - 1^ 
D-i 



(3.4) 
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For D = 5 and D = A, rg^cS = 2rs and rs,cff = (3\/3/2)rs, respectively. The surface area of the hole, Ajj, is 

Ad = nD-2ri-^, (3.5) 

where is the area of the {D — l)-diniensional unit sphere. The effective area Ad^^s is given by replacing the 

radius with rs^eff- 

Next we describe the Hawking evaporation process of braneworld primordial black holes. Recently numerical 
computations have been performed for emissions of matters into four dimensions in III and that of graviton into 
bulk 't^, "s^, "sT, '82^1 , which are extensions of the four-dimensional result by Page 0. For the emission of brane matter, 
the calculation is done with the 4D metric obtained by projection, where the deviation from the 4D result comes 
through the form of metric function. On the other hand, graviton emission is given using original 5D metric. The 
number of dimensions of momentum space is denoted as d, which we set 4 for brane matter or 5 for bulk graviton 
later. The differential emission rate for general particle species, indicated by j, is given by 

^^=^^ e^/^H±l(2.)^-i ' ^'-'^ 

where gj is the number of internal degrees of freedom of the particle, E — (fc^+m|)^/^ is the energy, aj is the absorption 
cross section depending on hole's mass, particle's mass, energy, and spin. The sign in the denominator expresses which 
statistics the particle obeys. Numerical results showed that in the high-energy regime, i.e., when the geometric optics 
approximation would be applied, the absorption cross section is given by aj ~ Ad^cff^d-3/[{d—2)Qd-2\, but in the low- 
energy regime a certain amount of deviation exists |77l [79l ISOl Isil Is^ . The total mass-loss rate, Stefan-Boltzmann's 
law for black hole, becomes 

dA4h\ _ S-Tp'^'^^dF- .(T^.A T'^- r3 7^ 

^Jevap -'\L " -lo{Tu>dA,^A (3.7) 

where ad is a d-dimensional analog of Stefan-Boltzmann constant given with gamma function T{d) and zeta function 
C(d) as 



and 7d(2h) is the effective degree of freedom fixed after numerical result, which is absorbed into e for convenience 

^_ (d-2)/(d-3) -^x (d-2)/2 /^_^\d 



e(TH) = 7.(Th)1^.-2.. j . (3.9) 

To illuminate the effect of dimensionality, we rewrite it as 

e= 1.118 X 10-^ Y.9dATii)fsCDA,s, (3.10) 

s 

where gd,s is the effective number of dof for each spin s depending on the temperature, /^'s are those introduced in 
[Tol | expressing (relative) spin-dependent deviation in four dimensions from black-body emission, and CD.d,s^s, are the 
ratio of emission power relative to that in 4D. We use Table II-III of [T^I for CD,d=4,s^2 to describe matter emission 
into brane and Table I of (SOl] for CD,d=5,s=2, bulk graviton emission^. Parameters we need are listed below 

= uncharged , = 0.060, = 0.007, 

10.142 with unit charge 

2 

CD=5.d=4:,s=l/2 — 14.2,C_D=5 d=4.s=i = 27.1, C£)=5.c;=5,s=2 = 7 ^ 103. (3.11) 

■ 5 

The numerical factor 1.118 x 10~^ keeps the four-dimensional normalization of fg's such that massless dof gives 
2 X fs^i X 1 4- 6 X fs=i/2 X 1 = 1, which actually recovers (dA/bh/dt)cvap = —5.34 x 10^^(Afbh/l g)~^ g s~^ [l3|. Since 



^ The values for graviton in 80] includes polarization dof. 
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we are now in the braneworld, the degrees of freedom on the brane, namely the standard model particles, takes 3D 
phase space of momentum and black hole surface area is 2D, while that in the bulk, graviton, takes 4D phase space 
and the area is 3D. Thus the mass- loss rate is formally given as 



(3.12) 



evap 



where g4 is calculated using the degree of freedom of the standard model matter, while 55 is the degree of freedom of 
graviton polarization reflecting dimensionality as ,9,5=2. d = d{d — 3)/2. The above formula is apparently the same as 
in four dimensions except the added five-dimensional term for graviton dof. However, the Schwarzschild radius does 
no longer obey four- dimensional relation for those holes much smaller than the compactification scale £. For a small 
enough black hole, 



dAi, 



bh 



dt 



evap 



Ma 



M4 



where the numerical factor is calculated according to the temperature of the hole as 



a5 = — e : 



0.0227 7, i/'s, i/'s, 5( 
0.0333 7, i/'s, i/'s, fif, e, e+ 



(3.13) 



(3.14) 



Through a period in which as is regarded as constant, the solution is obtained as 

1/2 

/ ;\//,_ = /ex ' I f — f, \ 



V M4 



2a. 



t-ti 



Mi 



(3.15) 



where Mbh.i = Mbh(ii) is introduced as some initial condition. The duration of this phase can be estimated as 

J~4 



At- 



1 

2^ 



V M4 



For comparison, we also show the 4D relations below 



/ dMbh \ 



= — Q4 



/Mbhy'M4 



a4 = 



Mbh(t) 
At 



V M4 



0.000283 7, !/'s, i/'s, 5 
0.000442 7, z/'s, i/'s, 5, e, e+ ' 

-,1/3 

/ t-t: \ 

3a4 



t-ti 



t4 



M4, 



1 

3a4 



M, 



bh,i 



M4 



(3.16) 

(3.17) 
(3.18) 

(3.19) 
(3.20) 



As we will discuss photons emitted from currently evaporating PBHs, our attention is concentrated on sufficiently 
large £ cases in which there exist 5D PBHs with lifetimes comparable to the present age of the Universe. For a typical 
5D black hole with lifetime to = 13.7 Gyr, we can assume their temperature is low enough so that its emission channel 
is restricted to massless species and, at most, e/e+. Then its primordial mass Af^j^ p is calculated as 



Mb*h,p 



Its corresponding temperature is 



2a5 



TjJ ~ 50 



1/2 



M4 ~ 5 X 10^ 



0.1 mm 



-1/2 



0.1 mm 



-1/4 



keV, 



(3.21) 



(3.22) 
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from which we confirm the consistency of low-energy approximation. The dominant contribution to the currently 
observable photon spectrum comes from this Af^j^ p holes. The condition for the five dimensionality of M^j^ p PBH is 
set as 

rs(A4Vp)^10-^(^^^^) ' mm<L (3.23) 

This gives £ > 10^^^ m. On the other hand, such 5D PBHs can be formed only when the reheating temperature is 
sufiiciently high such that at least /Mh i-h < -^^bhp- This condition reads 

/ / \-5/16 

Tr-h > lOV'/' TH GeV, (3.24) 

\U.i mm / 

which is, for large £, stronger than Eq. (|2.5U|I . The presence of accretion, F > 0, makes this lower limit stringent. 
F = 1 with £ = 0.1 mm and / = 1 gives T,h > 10^ GeV. In the 4D case, we get T,h > lO^/^/^ GeV for the same 
argument. 



B. Constraints from the extragalactic diffuse photon background 

It is not the purpose of this paper to show constraints on perturbation spectrum for the whole region of the 
parameter space. Rather, we discuss the most remarkable case, namely £ = 0.1 mm, and see how the existence of the 
extra dimension modifies constraints on the perturbation from four dimensions. In addition, two distinguishable cases 
for accretion efficiency, F = 1 and F = 0, will be examined. We give constraints on the set of the spectral index n 
and the reheating temperature Trh- We calculate the photon surface brightness for each set and determine if the set 
is allowed or not by the criterion that the whole spectrum is below the observed data or not. The photon background 
data used here are taken from [s^, 0j EE Ell (we also show the gamma-ray result in 87] ) . 

Let us formulate the diffuse photon spectrum from PBHs. At some fixed time t, the photon emission power is given 

as 

where (1 -I- z)^ expresses the decreasing rate of number density. For the emission rate per hole, we take = 2 and 
d = 4 for Eq. H3.6|l and use geometric optics approximation 

7 _ f,oH ^3 26) 



dEdt e^/^H _ 1 

Here we ignore secondary photons. Photons suffer red-shift and the energy that we will observe is Eq — E / {1 + z). 
Energy density U decreases as (1 + 2)"* to give its present density Uq 

An observable quantity which we use is surface brightness given as 

K^Eo) = JLE^iEA keV cm-2 sr"! keY-\ (3.28) 
An Eq 

We observe that there are two possibilities one should take care: exponential spectrum or power-law spectrum 
around M^^ according to the spectral index n > n* or n < n*, respectively, where n* is given by equating 5h,min ~ 
ah(M4) as 

log[<5h,mi„/(3/5)ao] _ 4.0 



log 



8.4-^3.0?/+ (0.014 + 0.12?/) log(£/0.1 mm) + (0.08 + 0.08y) log/' 



(3.29) 

where y = F/{w — F). In the case oi £ — 0.1 mm and / = 1, this gives n* « 1.5 for F = and n* « 1.4 for 
F = 1. First we examine the possibility of n > n*. In general such a large n causes excess above the observed photon 
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background, so it appears that n > n* is rejected. However, we should remember that the hghtest mass of existing 
PBHs, Mbh.min = /-^^h,rh, IS determined by reheating temperature. If Tih drops below ^ 100/^/^(^/0.1 mm)~^/^^ TeV, 
then the mass of the lightest PBH goes above M^j^ (see Eqs. (j2.49|l and (j^2J) and the constraint from the photon 
background makes no sense. On the other hand, n <^ n* case is not constrained any more since high reheating 
temperature only increases lighter PBHs than M^^i p which virtually does not contribute to the photon background. 
The problem reduces to a question how far the upper limit of n can approach n*. Some illustrative examples of 
the photon spectrum are shown in Fig. UTTbI The case with the large extra dimension, = 0.1 mm, with two 
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FIG. 2: Illustration of constraints on the photon spectrum. Each panel corresponds to a combination of £ and F. For each n, 
one of the spectra for lower Trh is allowed and the other for higher Trh is rejected. 



possibilities of the accretion efficiency, F = 1 and F = 0, are considered. For each set of {i,F), four photon spectra 
for different perturbation spectral index n are shown; two are allowed by either smaller spectral index or lower 
reheating temperature, while the other two are rejected for either larger index or higher temperature. From these 
figures, we observe that the constraint on the spectral index is the strongest for {£ — 0.1 mm, F = 1). This is because 
the mass spectrum is the hardest for this case so that the holes with mass around M^^^ are potentially most abundant. 

Fig. IIII Bl shows that every spectrum in the 5D cases makes a peak in the energy region of the High Energy 
Astrophysics Observatory (HEAP) J83l|. One notices that the spectra are significantly red-shifted compared to those 
presented in our previous work |37|, which were rather similar to black body emission with temperature TjJ. This 
refiects two facts; one is that the spectra in IST^ were drawn assuming scale-invariant power spectrum, whereas it is 
turned to be not allowed in the present setup from the above arguments. An exponentially decaying mass spectrum 
leads to a relative reduction of the population of heavier PBHs than M^j^. The other comes from lowered reheating 
temperature, which reduces the number of lighter PBHs than M{Jjj. Both the two facts indicate that the contribution to 
photon spectra from holes other than M^^ ^ is enhanced, hence the superposed photon spectra are seen to be relatively 
red-shifted. Note that even if the perturbation index n is more blue, say n ~ 2, lower reheating temperature could 
save the situation by virtue of this "reddening" effect. Of course this is only the consequence of the assumptions in our 
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method, namely, the single power-law density perturbation spectrum with the large-scale WMAP-|-SDSS-gal-f LyaF 
normalization. 

At last, the allowed region for the set (n,Ti.h) is given in Fig. lmBl As was mentioned previously, any n could be 
allowed no matter how it is large if only the reheating temperature is low enough. However, we here focus on the 
most meaningful region, that is, around the possible upper bound on the index n in the case of marginally high T^h. 
Approximately it is sufficient to require either conditions on Tih or n is satisfied. The constraints on the braneworld 



T 




- 1=0.1 mrn F=1 — -- ^ 1=0.1 mm F=0 

1.260 1.280 1.300 1.320 1.340 1.360 1.380 
Spectral Index n 

FIG. 3: Allowed regions of (n, Trh) for each combination of {£,F). Left-hand side is allowed and the opposite is rejected. The 
4D result is also shown in the inner panel. 



inflation for £ = 0.1 mm is summarized in Tab.^ Although only the case of £ = 0.1 mm has been considered, we can 

TABLE I: Requirements for the reheating temperature T^h and spectral index n. 
Dimensionality /Accretion Threshold To realize larger n: 
^ = 0.1mm/F = n « 1.35 Ta < 10^ GeV 

£ = 0.1mm/F = l n f» 1.28 T^h < 10^ GeV 

4D n f» 1.39 < 10* GeV 



anticipate what general £ results in. Considering Eq. H3.29|l . we deduce that both smaller £ and smaller F will weaken 
the constraint. 

We also obtained an updated constraint on four-dimensional PBHs. Since we used the observational data of the 
small-scale density perturbation obtained by WMAP, SDSS, and Lya Forest, the constraint on the spectral index n 
is relaxed compared to the result which used the COBE data normalized at the present horizon-scale. On the 
other hand, the border of the reheating temperature is unchanged. 



IV. DISCUSSION AND CONCLUSION 



In this paper, we obtained the primordial black hole mass spectrum produced from inflationary density perturbation 
in the Randall-Sundrum braneworld. The result is Eqs. (|2.43|l and 1)2. 47|) . Fig. ^ showed typical shapes of the mass 
spectrum. Apart from the normalization, the mass spectrum contains two scales. One corresponds to the scale of 
the transition from the p-square to the ordinary radiation-dominated era. Below this scale the mass spectrum is 
deformed. In the case of no accretion, F = 0, the index is softened to —17/8 from the four-dimensional value —5/2. 
On the other hand, in the presence of accretion, F ~ 1, the index is hardened to —2.65. The other scale in the mass 
spectrum is of the transition from power-law to exponential decreasing specified by ^h.min ~ crh(Mth). 

Using above results, we calculated the diffuse photon spectrum emitted from PBHs to constrain the perturbation 
spectrum n and/or the reheating temperature Trh- We considered a power-law power spectrum of curvature pertur- 
bation V-ji^ = A{k / kti)^~^ which is normalized at the Lyman-a forest scale fcn = 0.9/i Mpc~^ to be A = 2.28 x 10^^. 
Recent numerical results on the Hawking radiation in higher-dimensions |73, III) IZ^ were incorporated. 
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The observed photon background placed an upper limit on the photons from PBHs. Since the contribution to the 
photon spectrum is dominated by currently evaporating holes, the constraint was set so that the amount of such 
critical holes be sufficiently suppressed. Fig. If 11 Bl showed illustrations of spectra. By virtue of the rapidly varying 
form of the mass spectrum, the result presented in Fig. 1111 51 and Tab. were summarized simply: The requirement 
is either the spectral index n is smaller than some critical value, or the reheating temperature Trh is lower than some 
critical value. This upper limit for n could be estimated by n* in Eq. (|3.29|l . above which the mass spectrum around 
M^jj p becomes power-law, while the threshold for T^h was given from the condition /Afh(rrh) = -^^bh pi below which 
-^bh,p holes cannot be produced. It is interesting to note that the values of Tlh for five dimension cases in Tab. ^ 
is comparable to those values obtained for unstable gravitino decay |89l lOOl l9l] | , which are not connected with the 
dimensionality. 

It was helpful that each formula has been given a form which manifestly shows the recovery of the standard four- 
dimensional version ("p^" part) and the deviation caused by the higher dimension ("p^" part). Taking the recent 
observational data into account, we also have successfully updated the four-dimensional result on this subject obtained 
in the literature p^ . 

In the formulation, we owed a simplification with respect to the generation and evolution of super-horizon curvature 
fluctuation in the braneworld cosmology. There are still difhculties to handle the perturbation on the brane considering 
the bu lk ge ometry, but some intuition can be made from the results of tensor perturbations |93 . 193 . |9J, |9a, |9a, |93, 
IqM I99I llOC| . Gravitational waves produced in the inflationary regime or matter gravitational collapse can escape into 
the bulk, so the energy-loss into the bulk can also takes an important role in the formation process of primordial 
black holes. In this sense, omitting the perturbation of bulk geometry might lead us to an over estimate of the PBH 
production, and the constraint obtained in this paper may be weaken after the effects will be appropriately taken into 
account. 

The inflation mechanism itself in the braneworld scenario is rather unveiled. There are some suggestion for the 
brane inflation and generation of the densi ty perturbation on the brane |39l l4lj| , and we have had some constraints 
on the spectrum in the literature |lOlL llO j |. PBH offers an alternative way to investigate brane inflations. 
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